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Abstract 

The Dirac structure of the nucleon self-energy in symmetric nuclear matter as well as neutron 
matter is derived from a realistic meson exchange model for the nucleon-nucleon (NN) interaction. 
It is demonstrated that the effects of correlations on the effective NN interaction in the nuclear 
medium can be parameterized by means of an effective meson exchange. This analysis leads to a 
very intuitive interpretation of correlation effects and also provides an efficient parametrization of an 
effective interaction to be used in relativistic structure calculations for finite nuclei. 

1 Introduction 

Microscopic studies on the bulk properties of nuclear systems have shown that two ingredients may 
be needed to derive nuclear properties from realistic nucleon-nucleon interactions: the consideration of 
correlations beyond the mean-field approach and the relativistic structure of the nucleon self-energy in 
the nuclear medium. In fact it is known for many years that attempts to use realistic NN interactions, 
i.e. models for the NN interaction which have been fitted to the NN scattering phase shifts, in a simple 
mean field or Hartree-Fock calculation lead to nuclear systems which are unbound (see e.g. the recent 
review jy). Therefore, various techniques have been developed to account for correlations beyond the 
mean field approach including the Brueckner hole-line expansion^ |[, the coupled cluster or "exponential 
S" approach^, ||], the self-consistent evaluation of Greens functions ||, variational approaches using 
correlated basis functions j?], ||] and recent developments employing quantum Monte-Carlo techniques^, 

In the framework of the Brueckner theory the effects of two-nucleon correlations are taken into account by 
evaluating an effective interaction, the so-called G-matrix. This G-matrix corresponds to the T-matrix 
of NN scattering, however, in the nuclear medium accounting for Pauli and dispersion effects. It is 
obtained by solving the Bethe-Goldstone equation or in the case of relativistic meson exchange models 
for the NN interaction by solving an equation which corresponds to a three-dimensional reduction of the 
Bethe-Salpeter equation like the Blankenbecler-Sugar or the Thompson equation JIT|. In the Brueckner- 
Hartree-Fock (BHF) approximation the nucleon self-energy or single-particle potential is then evaluated 
in terms of this G-matrix. 

Accounting for the two-nucleon correlations in this way or by means of other many-body approaches, one 
obtains results for the saturation property of nuclear matter or the binding energy and radius of finite 
nuclei]!]], which are quite reasonable. All such results, however, form the so-called Coester band[]l2||, i.e. 
they either predict a binding energy which is to small or a saturation density which is to large (a to 
small radius in case of finite nuclei) as compared to the empirical values. Extensive studies have been 
performed to account for three- nucleon correlations |k| . One finds, however, that the inclusion of three- 
nucleon correlations yields a small effect only, the phenomenon of the Coester band persists for studies 
within the framework of Brueckner theory as well as within other approaches to solve the many-body 
problem. 

It has been suggested to consider three-nucleon forces, which are adjusted to produce the empirical 
saturation point of nuclear matter |Q. Another possibility to shift the calculated saturation point away 
from the Coester band towards the empirical value is to account for relativistic effects. These studies 
have been motivated by the phenomenology of the Walecka model[^5|. The NN interaction in this model 
is described in terms of the exchange of a scalar meson, er, and a vector meson uj. Calculating the nucleon 
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self-energy, S, from such a meson exchange model within a Hartree approximation, one finds that the 
u) exchange yields a component X , which transforms under a Lorentz transformation like the time-like 
component of a vector, while the scalar meson exchange yields a contribution X s , which transforms 
like a scalar. Inserting this self energy into the Dirac equation for a nucleon in the medium of nuclear 
matter leads to single-particle energies, which are as small as the empirical value of -50 MeV. This small 
binding effect, however, results from a strong cancellation between the repulsive X° and the attractive 
X s component. The attractive scalar component X s leads to Dirac spinors for the nuclcons in the nuclear 
medium, which contain a small component significantly enhanced as compared to the Dirac spinor of 
a free nucleon. This effect is often described in terms of an effective Dirac mass M * for the nucleon, 
which can be of the order of 600 MeV in nuclear matter around saturation density. This implies that 
the Dirac spinors for the nucleons in the nuclear medium are quite different from those in the vacuum. 
Since these Dirac spinors are used to evaluate the matrix elements for the meson exchange model of the 
NN interaction, this leads to a medium dependence of this interaction, an effect which influences the 
calculated saturation property. 

Also realistic models for the NN interaction contain large contributions from the exchange of scalar and 
vector mesons. For such meson exchange potentials, V, one can determine the Dirac structure of the 
nucleon self-energy using the Hartree or Hartree-Fock approximation in a straight forward way. As we 
have discussed above, using such realistic NN interactions, one has to account for correlations beyond 
the Hartree-Fock approximation and determine the self-energy in terms of the G-matrix rather than the 
bare interaction V, i.e. perform what is called a Dirac-Brueckner-Hartree-Fock (DBHF) or a relativistic 
BHF calculation [l?], [lj, [l9| ^0|. Since, however, this G-matrix is obtained as a solution from a non- 
relativistic reduction of the scattering equation, it provides matrix elements only between single-particle 
states and does not keep track of the relativistic structure of the effective interaction. 
Various approximation schemes have been developed to determine the Dirac structure of the self-energy 
X or the structure of the nucleon Dirac spinors within the context of the DBHF approach. A rather 
simple scheme has been suggested by Brockmann and Machleidt Jl7| . They determine the momentum 
dependence of the single-particle energy. Identifying this single- particle spectrum with a corresponding 
spectrum derived within the Dirac-Hartree approximation one can extract the effective Dirac mass M*. 
The underlying assumption is that the nucleon self-energy is dominated by the scalar, X s , and time- 
like vector component, X°, which are constants independent of the momentum of the nucleon. This 
approximation seems to work reasonably well for symmetric nuclear matter but it fails in the case of 
neutron matter pl|. Therefore calculations of the equation of state for asymmetric nuclear matter, which 
are based on this approach p^, p3{ , should be considered with some caution. 

Another scheme, the so-called projection method, analyses the antisymmetrized matrix elements of G in 
terms of sets of operators, which are invariant under Lorentz transformation. If the relativistic structure 
of the effective interaction G is defined in this way, one can derive the Dirac structure of the self-energy 
and determine the density dependence of the nucleon spinors |24[ |25| . As will be discussed more in detail 
below, the result of this analysis depends on the choice of relativistic invariants. As an example we 
mention the well known feature that the one-pion-exchange contribution yields identical matrix elements 
for the positive energy Dirac spinors using pseudo-scalar or pseudo- vector coupling, the Fock contributions 
to the scalar part of the self-energy, however, are quite different^, ^(| ^7j. As another example we will 
discuss the p exchange below. 

In order to minimise these uncertainties of the projection method, we suggest to split the G matrix into 
the Born contribution, the bare interaction V, and the corrections due to the correlations. While the 
Dirac structure of V is well defined one may employ the projection method for the correction term only. 
It turns out that these corrections can be described rather well in terms of the exchange of a few effective 
mesons with high masses, reflecting the short range of the NN correlation effects, and coupling constants 
depending on the nuclear density. The results of this analysis provides some insight into the effects of NN 
correlations. This analysis is similar to previous attempts by Boersma and Malflietj2i| and Elsenhans 
et al.|2^]. Because of its simple structure, however, it might be more appropriate to be used in DBHF 
studies of finite nuclei p9|, |3(| 

After this introduction a brief description of the main ingredients of the projection technique will be 
presented in section 2. Results of the analysis for the effective NN interaction in symmetric nuclear 
matter and neutron matter derived from various models of the Bonn potentialjll]]. After a detailed 
discussion of the results in section 3, the main conclusions will be summarised in section 4. 
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2 Relativistic structure of the G-matrix 



The Dirac equation for a nucleon with momentum k in a medium of nuclear matter can be written 

{jk — M — T,(k))u(k, s) = (1) 

with the self-energy E(fc) accounting for the mean field generated by the nuclear medium. By the 
requirement of translational and rotational invariance, parity conservation and time reversal invariance, 
the general form of the Dirac structure of the self-energy is given in the nuclear matter rest frame as 
4 x 4-Matrix by 

E(jfc) = £ s (fc) -7°£°(fc) +7- fcE"(fc) (2) 

where £ s , E°, £" are functions, depending for on-shcll nucleons (k° — E(k)) only on the absolute value 
of the three- momentum k = \k\ and the Fermi- momentum kf, which is related to the density via p = 
S/(3iT 2 )k 3 f where the isospin degeneracy yields 6 = 2 for nuclear matter and 6=1 for neutron matter. 
The density dependence will be suppressed throughout this section. 
The components of the self-energy are easily determined by taking traces of the form 



E s (fc) = iTr[E(fc)], E°(fc) = -lTr[ 7 °E(fc)] , E«(fc) = --TVtf • fc£(*)] C!) 



Introducing the effective quantities 



jvr(fc) = M + ^ (fc) E *{k) = ^M (4) 

1 + E" (k) v ' 1+s W 



the Dirac equation can be rewritten in compact form setting k* = k 

(jk* - M*)u(k, s) = (5) 

which is formal identical to the Dirac equation in the vacuum case. Therefore the positive energy solution 
to eq. (||) is given as 

,? *, lE*(k) + M*(k) ( 1 \ 



with the covariant normalisation 

' " " _ M*(fc) 



u« = 1 , u+m = (7) 



and the in medium on-shcll relation 

£T 2 (fc) = M* 2 {k) + k 2 (8) 

The self-energy itself is connected to the two-particle effective interaction in the medium, the G-matrix, 
through the Hartree-Fock relation 

f 

Ea/j(fc) = —if J7r-T79ra{q) {G(k,q;k,q) acr . f j T ~ G(k,q;q,k) aa . Tl3 } 

Jq<k } y^r 



d 4 q . 



g T a{q)G (k, q; k, q) aa -p T 

q<k f 



"' q Ti-[g(q)G A (k,q;k,q)] a p (9) 



q<k f 



(27T) 4 



with the propagator for real nucleons propagating on-shell inside the Fermi-sea in the nuclear matter rest 
frame 

9rAq) = + M*) T J(q° E(q)Mk f - \q]) (10) 

The equation above (^) must be understood as an operator equation. Therefore, also the G-matrix is 
needed in terms of Dirac operators as will be discussed now. 

One of the main ingredients of the DBHF approach is the derivation of the effective interaction G in 
the medium from a realistic nucleon-nucleon interaction V. This is achieved by solving the relativistic 
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Bethe-Goldstone equation, a kind of Thomson equation for the scattering of two nucleons in nuclear 
matter 

- - F tfik - -. M* 2 (l(k P\ 

G(q>,q\P,u>)=V(q>,q) + / —-V{q>,k) ^>> G(k,$P,u) (11) 

J ^ ' P+k P+k 

with the starting energy 

- = 2Et +3 (12) 

where P = h(pi +P2) denotes the center of mass momentum and q = \{j>\ ~ P2) is the relative momentum 
of the initial state, k and q' are the relative momenta of the intermediate and final states of two interacting 
nucleons. The accessible intermediate states are restricted by the Pauli operator Q to momenta above the 
Fermi momentum kf . To shorten the notation isospin indices will be suppressed throughout this section. 
To account for the spin, the G-matrix (A^A^G^', (?) | Ai A2 ) is conveniently obtained in a basis of helicity 
states, where A^ , X' 2 and Ai , A2 are the helicities of the final and initial state of particle 1 and 2. Demanding 
parity conservation, time reversal invariance and conservation of total spin reduces the 16 possible helicity 
amplitudes to six independent amplitudes, which are normally chosen to be 

Gi = {++ \G(q>, q)\++), G 2 = <+ + \G(q>, q)\ - -) , 

G 3 = (+- \G(q',q)\ + -) , G 4 = <+ - \G(q',q)\ - +) , 

G 5 = {+ + \G(q>,q)\ + -), G 6 = (+ - \G(q', q)\ + +) , (13) 

where +, — denotes the helicity | and —5, respectively. The explicit dependence on P and ui has been 
omitted. The number of independent helicity matrix elements is further reduced to five in the case of 
on-shell scattering, \q'\ = |<f|, which implies 

G 5 - -G 6 (14) 

In order to evaluate the nucleon self-energy eq. (||) , the operator structure of the G-Matrix is needed for 
the on-shell case in the nuclear matter rest frame; on the other hand the determination of the operator 
structure is most easily done in the center of mass frame where P — and the scattering angle D between 
the relative momenta q', q is fixed at 1? = 0. Therefore, only q = \q\ have to be considered in the 
following as argument. The operator structure is then given by expanding the G-matrix in terms of five 
independent Fermi covariants according to the five helicity matrix elements of the G-matrix: 



{>! 1 ,x! 2 \G(q)\x 1 ,X2)A = £ r^) ( |ff >ff l^-r^HiTfrf' 

i=S,V,T,A,PS 



X 



£ r A (q) (AUUl^Tflg^A^ (15) 



i=S,V,T,A,PS 



^{I^^tV,! 5 } (16) 



with the Fermi covariants 
and the Lorentz invariant amplitudes 

V^q) ^{q) -Y.Fn^xiq) (17) 

fc=i 

where the F^i is the well known Ficrz- Transformation 
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(18) 



The subscript A in eq.(|15|) indicates, that only antisymmetrized matrix elements are obtained as solution 
of the Bethe-Goldstone equation and therefore only those can be analysed. Here, the explicit splitting of 
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the antisymmctrized matrix element in its direct and exchange part (labelled with the subscripts D and 
X) only illustrates, that, using the local Fermi covariants eq.(|l6|), both parts can be rewritten with the 
help of the Fierz-Transformation in terms of the antisymmetrized amplitudes Y\ and the direct matrix 
element of the Fermi-covariants. Nevertheless, the direct and exchange part cannot be determined from 
the antisymmetrized matrix elements and it is not required to evaluate the self-energy. 
Inverting eq.([l5]) yields the antisymmetrized amplitudes TV where only the scalar and vector T\ 
amplitudes are needed to evaluate the self-energy using eq(H). 



^(fci) = -^j o f d 3 k 2 ^T\{q) (19) 

Here, k\ and k 2 are the single particle momenta in the nuclear matter rest frame. The relative momentum 
q, defined in the cm. frame, is related with fci and k 2 via q = y^s/A — M* 2 where s — (E* (k\)+ E* (k 2 )) 2 — 
(fci + k 2 ) 2 is the invariant mass. This transformation between the nuclear matter rest frame and the center 
of mass frame is described in detail inp4||. 

The used set of Fermi-covariants is sufficient to reproduce the matrix elements of the effective interaction 
G for the positive energy solutions of the Dirac equations. From the resulting amplitudes F A (q) and 
r^(g) one can immediately determine the components of the self-energy according to ([H]). Nevertheless 
this procedure is not unique and depends on the chosen set of operators T. A well known example which 
illustrates this dependence is the pion-exchange part of the NN interaction, which either is described by 
a pseudo-scalar or a pseudo- vector coupling. Both couplings yields for on-shell nucleons the same matrix 
elements, if the coupling constants f pv and g ps obey the relation f pv /m n — g ps /(2M), where m„ is the 
pion-mass and M the mass of the nucleon. Nevertheless, the components of the self-energy are complete 
different for both couplings. A more detailed discussion of this and other examples will be given in the 
next section. 

To circumvent aforementioned problem, which already occurs in the case if the effective interaction G is 
replaced with the bare nucleon-nucleon interaction V, the G-matrix is split into two parts 

G = V + AG (20) 

The decomposition eq.(|l5|) is applied only to the residual part AG since the explicit Dirac structure of 
V and the way to evaluate the components of the self-energy from V is known explicitly. 



3 Results and discussion 

In the first part of this section we would like to show some ambiguities of the projection method which is 
related to the choice of the covariant operators in (|l6|). For that purpose we consider various components 
of the bare NN interaction V using as an example the OBE potential Bonn B as defined in table A. 2 of 
[fj"l"l . If one considers only the 7r-exchange contribution and analyses the matrix elements of by means of 
the projection formalism discussed in section 2, one obtains large components for Lorentz invariant scalar 
and vector amplitudes T B A and T A which are due to the Fock-exchange terms of the pseudo-scalar operator 
representing the it exchange in the projection formalism. Integrating these amplitudes as indicated in 
( |l9[) leads to large scalar (£ s ) and time-like vector components 

This is visualised in the left part of Fig. [j], where these components calculated for symmetric nuclear 
matter at the empirical saturation density (Fermi momentum kp = 1.36 fm -1 ) are presented as a function 
of the nucleon momentum k±. Note that this analysis yields an absolute value for the scalar component 
which gets as large as 600 MeV, which implies that the effective Dirac mass, M* , of the nucleon in the 
nuclear medium would become as small as 340 MeV. Furthermore this analysis yields quite a strong 
momentum dependence of the scalar and vector components of the self-energy E. Because of the strong 
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cancellation between the scalar and vector components, the total effect of the 7r-exchange on the single- 
particle energy 



is rather small. 

Calculating the self-energy directly from the 7r-exchange part in the Bonn B potential leads to quite 
different results as one can see from the right-hand side of Fig. |l|. In this case the scalar part of the self- 
energy has the opposite sign, it is small (of the order of 15 MeV) and exhibits a very weak momentum 
dependence. The single-particle energy for the nucleons in nuclear matter, calculated according to (pi]), is 
identical to the one derived from the projection method, the Dirac structure of the self-energy, however, 
is completely different. The reason for this difference is well known: The matrix elements of the 7r- 
exchange potential, calculated for the positive energy solution of the Dirac equation for the nucleon, 
are independent of assuming either pseudo-scalar or pseudo-vector coupling for the 7rN vertex. The 
coupling between the positive and negative energy spinors, however, is quite different. The projection 
method analyses the matrix elements of the positive energy solutions in term of a pseudo-scalar operator, 
leading to large components of the self-energy. Assuming a pseudo-vector coupling for the pion, as it is 
done in the original NN interaction V yields much smaller components in the self-energy. This feature 
is known for a long time (see e.g. |Tq| ). It is the origin of the strong momentum dependence of the 
self-energy observed in pq , where the projection method has been applied. The problem might be 
cured by analysing the interaction in terms of the pseudo-scalar operator, transforming the pseudo-scalar 
component into a corresponding pseudo- vector term and evaluate the self-energy with this pseudo- vector 
part (see e.g.||j| |27||). 

Such ambiguities of the projection method, however, not only arise for the 7r-exchange part. Another 
example is of course the exchange of the pseudo-scalar - isoscalar meson, the r\ meson. Because of 
its weaker coupling constant and its higher mass, it does not play such a significant role as the tt- 
exchange. Nevertheless, also in this case the differences between the self-energy components derived via 
the projection method and from the direct evaluation are non-negligible as one can see from Fig. El 
The /5-exchange contribution to V is considered as a last example for such a comparison between the 
self-energies calculated via the projection method and the direct evaluation of the meson exchange term. 
Results are displayed in Fig. |^. The original /9-exchange term of the Bonn B potential contains a vector 
coupling but also a strong tensor coupling (Pauli-coupling) term in the Lagrangian for the pN interac- 
tion. Also these two different coupling modes cannot be resolved by means of the projection method. 
The difference between the direct evaluation of the self-energy and using the projection formalism are 
remarkable also in this case. 

Therefore, in order to minimise the ambiguities in the projection method, we suggest to split the effective 
interaction G according to eq.(|20|) into the bare interaction V and the correction term AG, representing 
the corrections which are due to the correlations. The Dirac structure of the bare interaction is directly 
known, and the projection method has to be applied to the analysis of AG only. This scheme is advan- 
tageous also from another point of view: One must keep in mind that the matrix elements of V and G in 
the helicity basis are dominated by the one-pion-exchange contribution. Note that the ratio of coupling 
constant and mass of the meson, g^/rn 2 al which is a measure for the importance of the various meson 
exchange contribution is about a factor 25 larger for a = it than for the a or u> meson. Analysing AG 
this dominating 7r-exchange contribution, including the momentum dependencies, which are related to 
the form-factors for the 7rN vertex, are removed, which stabilises the numerical analysis significantly. 
Results for the nucleon self-energy calculated again for symmetric nuclear matter at a Fermi momentum 
kp = 1.36 fm _1 are displayed in Fig. |[ The contribution of the bare NN interaction V to the scalar and 
vector part of the self-energy (see dashed lines in both parts of the figure) yields rather similar values 
ranging from -420 MeV to -400 MeV for k\jkp between and 1. This implies that the single-particle 
energies tend to be positive and one obtains no binding energy for nuclear matter. This is in line with 
the non-relativistic studies mentioned in the introduction: Using realistic NN interactions one does not 
obtain any binding energy if correlations beyond the Hartree-Fock approach are ignored. 
Adding the contributions of AG to the various components of the self-energy, one obtains the results 
displayed by the solid lines in Fig. [|. The correlation effects contained in AG reduce the absolute value 
for the scalar as well as the vector component of the self-energy. One could argue that the short-range 
correlations lead to a reduction of the wave function at small relative distances. This quenching of the 
relative wave function reduces the effects of the a exchange, which is the dominating contribution to E s , 
as well as the oj exchange, which is the driving term in E°. Since the ui exchange is of shorter range 




(21) 



G 



than the a exchange, the mass of the u) meson (m w = 783 MeV) is larger than the one of the a (m CT — 
550 MeV) , the quenching of the short-range components in the relative wave function is more important 
for the ui than for the a exchange. This explains that the absolute value of E° is reduced stronger by 
AG than the absolute value of S s . This difference on the effect of correlations on E s and S° leads to 
attractive single-particle energies and binding energy of nuclear matter. 

Its worth noting that the inclusion of AG leads to a constant shift in the components of the self-energy, 
which is almost independent of the momentum of the nucleon. This very weak momentum dependence 
of the AG effects can also bee seen from Fig. 0, where the antisymmetrized amplitudes (see eq.(|l5|)) 
arc displayed. These amplitudes are derived from the analysis of AG. Also these interaction amplitudes 
are almost independent of the relative momentum q. Deviations from the constant value are observed 
only for q — > kp, i.e. for momenta for which the energy denominator in the Bethe-Goldstone equation 
approaches the pole. 

This observation suggests to parameterize the component AG in terms of the exchange of effective mesons 
with infinite mass. To put it in different words: The correlation effects contained in AG are described 
in terms of an effective interaction with zero range. If we focus the attention on the scalar and vector 
interaction amplitudes only, at each density 4 coupling constants are required to parameterize AG in the 
form 

AG = [g s . s \\ + g v ^l^ + b s , H + v/Tjr (22) 
The results for these effective coupling constants are displayed in Fig. g. 

The value of these coupling constants indicate that AG is in general weaker than V but yields a non- 
negligible correction to the bare interaction V. The parameters are rather similar for the three versions 
of the Bonn potential (Bonn A, B and C as defined in table A. 2 of (ll)). The absolute values are typically 
larger for Bonn C than for the other two reflecting the fact that this potential yields slightly stronger 
correlation effects. The density dependence of these parameters is weak but non-negligible. The largest 
absolute values are observed for vector - isoscalar pseudo meson (g v ,s)- This reflects the fact, which we 
mentioned already above, that correlations yield a suppression of the w meson exchange in particular. 
Results for the properties of symmetric nuclear matter are displayed in Figs. and ||. Fig. ^ displays 
the density dependence of the effective Dirac mass M* (k) defined in eq. (§). The shaded area at each 
Fermi momentum kp indicates the range of values for M*(k), which is obtained at the corresponding 
density using the projection method. For a comparison we also present the effective mass derived from 
the momentum dependence of the single-particle energy as proposed by Brockmann and Machleidt Jl7| . 
The analysis of the single-particle energy, which is much simpler than the projection scheme, yields rather 
similar results for symmetric nuclear matter. Therefore also the calculated binding energies, derived from 
the projection scheme (see Fig. ||) are very close to those obtained in p7[ . 

Significant differences between these two methods, however, can be observed in the case of isospin asym- 
metric matter, like e.g. pure neutron matter. The effective masses obtained from the single-particle 
spectrum are significantly larger than the range of M*(k), which are derived at the same density from 
the projection method (see Fig. |J). This confirms the observation of Ulrych et al. |^l[ who studied a 
parametrization of G published by Boersma and Malfliet (2^] . This demonstrates that the analysis of the 
the single-particle energy is not a reliable tool to determine the Dirac structure of the self-energy: acci- 
dentally it works well for symmetric nuclear matter for the NN interactions considered, it fails, however, 
for asymmetric nuclear matter. 

How does this affect the calculated binding energies. As an example results for the binding energy of 
pure neutron matter calculated for the Bonn A potential are presented in Fig. The two methods yield 
almost identical results at low densities. At densities around 0.4 fm -3 , which corresponds 2.5 times the 
saturation density of nuclear matter, the evaluation based on the simple analysis of e{k) underestimates 
the energy by 10 percent. 

Finally, we add a remark on the asymmetry dependence of the parametrization of AG according to (p2|). 
If one employs the density dependent parameters for AG derived in symmetric nuclear as displayed in 
Fig. |^, assuming that the parametrization of AG is independent of the isospin asymmetry, also for pure 
neutron matter one obtains effective masses (M* at k = kp) and contributions of AG to the energy per 
nucleon in neutron matter as presented by the dashed lines in Fig. These results are very close to the 
corresponding values derived from a direct determination of AG in neutron matter. This indicates that 
the underlying assumption, the parametrization of AG is independent of the isospin asymmetry is quite 
reasonable. Therefore one may use the parametrization of AG as displayed in Fig. |^ for Dirac-Brueckner- 
Hartree-Fock studies of finite nuclei with N = Z as well as N ^ Z. 
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4 Conclusions 

A method is presented which determines the Dirac structure of the Brucckner G-matrix from its matrix 
elements between positive energy spinors only. The usual projection method is model dependent as it 
depends on the choice of the covariant operators. This model dependence is demonstrated for various 
meson-exchange terms. In order to minimise the model-dependence of the projection method G is split 
into the bare interaction V and the correction AG reflecting the effects of NN correlations. Since the 
Dirac structure of V is known, the projection methods is applied to the correction AG only. This analysis 
allows an explicit study of the correlation effects on the Dirac structure of the nucleon self-energy. A 
simple parametrization of AG and its density dependence is presented in terms of the exchange of pseudo 
mesons with infinite mass. This parametrization could be useful for the study of finite nuclei. It is shown 
that the simple method, which determines the Dirac structure of the self-energy from the momentum 
dependence of the single-particle energy yields fairly good results for symmetric nuclear matter but fails 
for asymmetric matter. 

This investigation has been supported by the SFB 382 of the "Deutsche Forschungsgemeinschaft" . 
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Projection 



Exact 




Figure 1: Contribution of the it exchange to the scalar (S s ) and vector component (£°) of the nucleon 
self-energy in nuclear matter with kp — 1.36 fm -1 . Results derived from the projection method (pseudo- 
scalar coupling) are displayed in the left part of the figure, while those derived from the original Bonn 
interaction (pseudo- vector coupling) are given in the right part. Note the different scales in the two parts 
of the figure. 
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Figure 2: Contribution of the r\ exchange to the scalar and vector component of the nucleon self-energy 
in nuclear matter. Further details see figure |l|. 
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Figure 3: Contribution of the p exchange to the scalar and vector component of the nucleon self-energy 
in nuclear matter. Further details see figure [l]. 
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Figure 4: The scalar (S s , left part) and time- like vector component (E°, right part) of the nucleon 
self-energy in nuclear matter at kp = 1-36 fin -1 . Results obtained for the bare NN interaction V (Bonn 
B potential of dashed lines) are compared to those derived from the G-matrix (solid lines). 
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Figure 5: Antisymmetrized interaction amplitudes T l A , derived from the analysis of AG presented as a 
function of the relative momentum of the interacting nucleons in their center of mass frame 
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Figure 6: Coupling constants for the parametrization of AG according to eq.(p2|) for the three versions 
of the Bonn potentials. These coupling constants are presented as a function of the Fermi momentum 
fcp of symmetric nuclear matter. The values of the constants are multiplied with the square of the mass 
of the a meson (m a = 550 MeV) and divided by Air to obtain dimensionless quantities of the order 1. 
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Figure 7: Effective mass M* characterising the Dirac spinor of the nucleoli in nuclear matter as a function 
of the Fermi momentum kF- The range of the momentum-dependent masses derived from the projection 
method are indicated by the shaded area, the result obtained from an analysis of the single-particle 
spectrum is displayed by the solid line. The data have been derived from the Bonn A potential. 
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Figure 8: Calculated binding energy of nuclear matter as a function of the Fermi momentum kF using 
the three different versions of the Bonn potential 
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Figure 9: Effective mass M* characterising the Dirac spinor of the nucleon in neutron matter. Further 
details see Fig. [?] 
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Figure 10: Calculated binding energy of neutron matter as a function of the density using the Bonn 
A potential. The results derived from the projection method are compared to those derived from the 
analysis of e(fc). 
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Figure 11: Contribution of AG to the energy of neutron matter and the effective mass M* , calculated 
as a function of hp. The results of the direct calculation (solid lines) are compared to 
the prediction of the parametrization of AG according to eq.(|2^) assuming that this parametrization is 
independent of the isospin asymmetry. 
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